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ABSTRACT 

This  work  describes  a  method  of  determining  the  ratio  of  heat  and 
momentum  dif f usivities  from  micrometeorological  data  obtained  by  the 
Great  Plains  Turbulence  Field  Program  [/]  • 

The  method  involves  obtaining  first  the  value  of  the  Deacon  wind 
profile  parameter,  used  as  the  indicator  of  stability,  and  with  this 
value,  the  ratio  of  two  forms  of  the  Richardson  number,  which  is  theo- 
retically equivalent  to  the  ratio  of  heat  and  momentum  dif f usivities . 

This  investigation  was  carried  out  at  the  U.  S.  Naval  Postgrauuate 
School,  Monterey,  California,  during  the  period  February  1959  -  May 
1959,  in  partial  fulfillment  of  the  requirements  for  the  degree  of 
Master  of  Science  in  Meteorology. 

The  writer  wishes  to  express  his  appreciation  for  the  assistance 
nd  guidance  given  him  by  Professor  Frank  L.  Martin  in  this  L.     ga- 
tion. 
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1.   Introduction 

A  theory  is  introduced,  after  Lettau  [61,  and  Martin  Q33, 
which  permits  computation  of  special  forms  of  the  Richardson  number 
and  flux  Richardson  number.   These  equations  are  dependent  upon 
the  Deacon  wind  profile  parameter  «  ,  and  the  roughness  para- 
meter zQ.      Equations  for  the  computation  of  these  parameters 
are  derived  also. 

A  recent  theory  is  included  in  these  determinations,  and  in 
a  sense,  this  paper  constitutes  an  attempt  at  testing  this 
theory  and  comparing  its  results  to  that  of  the  many  previous 
investigators  in  the  field  of  atmospheric  turbulence. 

The  data  utilized  in  this  study  is  evaluated  and  is  considered 
to  overcome  many  of  the  considerable  limitations  inherent  in  data 
available  for  earlier  research. 

Correlation  coefficients  of  the  resulting  ratio  of  heat  and 
momentum  dif f usivities  with  the  two  forms  of  the  Richardson 
number  and  the  Deacon  wind  profile  parameter  were  computed^and 
the  dependence  of  the  ratio  upon  the  other  variables  is  discussed 
and  compared  in  the  light  of  knowledge  gained  by  earlier  investi- 
gators . 


2 .   Theory  and  derivations 

The  work  of  von  Karman,  Fraud tl  and  Nikuradse  has  she  . 
that  the  fundamental  relationship  for  the  wind  shear  within  a 
neutral  surface  layer  is  given  by  C2] , 

where  u  =  the  mean  speed  at  a  height  z  above  the  surface 

u^   =  y/^°/p    called  the  "friction  velocity"  in  the  neutral 
a         layer 

La_  s  the  shearing  stress  at  height  z 

p   =  fluid  density 

k   =  von  Karman 's  constant 

On  integration,  Eq .  (1)  gives  the  logarithmic  lav/  for  the 

neutral  surface  layer  (defined  as  the  layer  in  which  the  shearing 

stress  is  effectively  constant), 

which  applies  to  fully  rough  flow,  and  z    ,  the  roughness  parameter 
is  the  measure  of  surface  roughness. 

In  studying  the  wind  profiles  versus  the  logarithm  of 
height  as  ordinate  under  non-neutral  conditions,  Deacon  L2~\ 
noted  that  the  wind  profile  differed  from  the  logarithmic  profile 
in  that  it  showed  convex  curvature  under  unstable  conditions 
and  reverse  curvature  under  stable  conditions.   He  has  suggested 
that  the  relationship, 

should  exist  under  non-neutral  conditions,  where  A  is  the  wind 
profile  curvature  parameter. 


. 
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Applying  the  condition  that  u  =  0  at  z  =  z  ,  which  holds 
for  the  logarithmic  profile,  and  integrating,  this  expression 
becomes , 


L\. 


--  *ftKT- '}  C4-) 


Assuming  that  -rH*'  ..  remains  constant  in  the  layer  Zi  to 
zo ,  one  may  write  (2), 

Jit  -  rCVa.V"P-  O  ^ 

Determination  of  zQ  using  an  atmosphere  of  neutral  stability 
allows  one  to  calculate  (J  knowing  the  winds  at  the  heights  1 
and  2.   Eq .  (3),  with  i£i  regarded  as  known,  may  be  solved  for  6 
by  a  method  of  successive  approximations  known  as  the  Newton- 
Raphson  method.   In  this  method,  successive  approximations  to  a 
root  of  any  equation  f(y)=0  may  be  written  in  iterative  form  as, 


»*  =  m  -  m  (0 


f'Cf 

Rearranging  Eq .  (5)  into  the  form  f(y)=0,  where  y=l  -  6  ,  one 
obtains,  ,.  j. 

Placing  the  above  expression  into  the  iterative  form,  and  making 
use  of  the  identity,   j^&    I'M   -  lej. , 0  (2l. j  JLt^i O 

gives,  HW-SiV^  _      Mi     _      /'jg>i*ti  I  1 

U        -    ,■  U.UJ  tl.  U)     *'  J  A5\ 

*•"  " p    [ftfeJH&J  -  (^^(f)]  %<» 

This  method  of  obtaining  G  was  used  in  this  paper  and 
requires  only  several  iterations  to  obtain  the  accuracy  specified 
in  Section  4. 


Since  the  solution  for  (3  from  Eq .  (8)  is  dependent 
upon  the  roughness  parameter,  a  method  of  determining  its 
value  under  the  applicable  conditions  must  be  devised.   As 
indicated,  zQ   values  may  be  determined  under  neutral  condi- 
tions at  sunrise  and  sunset.   Since  z   is  indep 
lapse  rate   [3}  ,  use  is  made  of  the  following  dimensionless 
ratio,  ciue  to  Davidson  and  Barad  tl]  , 

R     -    ll^-a,)  +^v  iQJi    „   , 

into  which  the  expression  for  the  logarithmic  profile  may  be 
substituted.   This  gives 


ZUJ 


ZLt  inii +Jln2i  -k-lni-^  +in2,  -^Jnio 


(to) 


Upon  further  reduction  of  Eq .  (10)  an  expression  for  ;::..  is 


obtained , 


"2*  z   A^ 


-! 


2"JL  it  A  X  ^ 


>*. 


2fU 


(fO 


where    ^J^2Z  -  jm^  +  U^+l^i^Jhi  t  ^  &  ]^£  -UifJLi}  ~JU?^,i, 

a.   Turbulence  theory 

The  derivation  of  the  remainder  of  the  theory  in  this 
paper  is  dependent  upon  Lettau's  turbulence  model  [C)    , 
which  was  adapted  to  the  data  in  a  manner  similar  to  that 
of  Martin  f8]  . 

In  analogy  to  molecular  diffusion,  the  theory  of  turbu- 
lent mixing  defines  the  exchange  coefficient  as, 

"  (1  2) 


KjJ 


which  is  written  by  Lettau   C63  as 

4 


6  3) 


where  w  i   :a      Che  mixing  velocity  and  1  the  root -mean-square 
mixing  length  introduced  by  frandtl  (.3J .      Since  turbulence 
theory  also  states  that, 

t---fu™'=    A^Vc-,t)  (\H) 

one  may  write  that, 

X-    ?  ^^  05") 

where  u  is  called  the  virtual  friction  velocity. 

Combining  the  last  three  equations  we  obtain, 

Lettau  [oj   assume:;  under  adiabatic  conditions,  that 


Under  non-adiabatic  conditions,  Lettau  [0>J     makes  the  following 

(\8) 


assumptions,  <j  ,.■ 

Ot'"  -         WOr 


sre  Nfl=  Kii/K?4,  the  ratio  of  the  coefficient  of  eddy  diffusi-* 
vity  for  heat  to  that  of  eddy  diffusivity  for  momentum. 

Substituting  the  expression  for  w  from  Eq .  (18)  into 
Eq .  (19),  gives  upon  rearrangement, 

■A.         va^  ^^Yn    °  *■ 


It  should  be  noted  at  this  point,  that  in  the  last  term  of 
Eq.  (20),  the.  term  Ng  represents  the  efficiency  in  the  conversion 
of  thermal  energy  into  kinetic  energy.   The  last  term  of  Eq .  (20), 
which  will  be  denoted  x,  may  be  written 

where  X  is  a  particular  form  of  the  Richardson  number  introduced 
by  Lettau   f6^ .   Actually, 


X=R[($Ll  (2  2) 


but  X  will  henceforth  be  called  the  "Richardson  Number",  whereas 
x  will  be  called  the  flux  Richardson  number  C6^. 

Combining  Eqs .  (18),  and  (19),  Lettau   Oil  obtains, 

Combining  Eqs.  (3),  (15),  (lo),  and  (23),  enables  one  to 
derive  the  following  results, 

".--f^ffi)1"        <"> 

This  derivation  is  shown  by  Martin  \jj~}. 
(b)   Use  of  finite-difference  derivatives 

In  computations  with  Eq .  (25),  it  is  necessary  to  use 
finite-difference  derivatives..   The  instruments  were  arranged 
at  constant  logarithmic  increments  so  that  the  following 


transformations  due  to  Lettau   D^  are  of  value, 


where  2]0  represents  the  geometric  mean  of  any  layer  :.»  to  Z2  > 
and  n*Z,  for  the  data  levels  utilized  in  this  work.   Therefore, 

and  in  terms  of  finite-differences,   /$%-  becomes  , 

Letting  z=2j^   makes  Eq.  (28)  representative  at  a  mean  height 
in  the  layer  l,  to  Zo  >  and  ^*/^V  becomes , 


Rearranging  leads  to 


•  «*--&&  te*)"1'  fee) 

and  solving  for  x,  gives, 

x       -\\     £r^l*2   1  fe,.^      |1         (3i) 

Hence,  it  becomes  possible  to  solve  for  x,  knowing  w  ,   6,  and  z0, 
(c)  An  expression  for  wa 

Strictly  speaking  wa  varies  with  stability  and  with  geo- 
strophic  wind.  The  variation  with  stability  is  not  yet  accu- 
rately known  but  the  variation  from  one  adiabatic  time  to 
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another  may  be  assumed  to  be  linear.   An  expression  for  w   in 
terms  of  the  neutral- wind  profile  may  be  obtained  by  logarithmic 
differentiation  of  Eq .  (2),  with  ui<a   replaced  by  wa ,  giving, 

whicl , on  expansion,  becomes 

J^  z  ^-u,,^^]  =   m,    £33) 

where  /\  Ljp^  is  notation  for 

zA      zq 

and  In  vr  +  ir  — —  =  Aln2  with  instrumentation  as  set  up  in  the 
Al      Z2 

Great  Plains  project. 

Substitution  of  Eq .  (33)  into  Eq .  (31  gives 

Variation  of  Lu  i    linearly  to  coincide  with  the  bi-hourly  observa- 
tions between' the  neutral  times  at  sunrise  and  sunset  will  give 
representative  values  of  the  flux  Richardson  number  x  during 
the  day. 
(d)  A  convenient  formula  for  X 

The  particular  form  of  the  Richardson  number  used  in  the 

Lettau  theory  was  shown  in  Eq .  (21).   This  form,  henceforth 
g 

denoted  X,  may  be  written  as   C73, 


*=-§&  C35) 


Similarly  to  Eqs .  (26),  we  have  the  transformation 


■  *'=*&--  4fc  (3  6) 


. 


and  therefore  in  terms  of  finite-differences,  the  following  rela- 
tionships  hold  at  3— £i„-z.« 

Substituting  Eqs.  (37)  into  Eq .  (35)  gives, 

where  ^ao  -  Ua*    may  be  averaged  over  all  successive  pairs  of 
levels  to  be  A**^  A  of  Eq .  (33).   In  Eq .  (38)  9m  may  be  approxi- 
mated by  0lu  =  T2 .   Making  the  substitution  U-.2    ~   ^a.>  =  ^#aA  gives, 

Since  "   no  longer  appears  in  either  x  or  X,  but  bas  been 
replaced  by   AUa,  the  latter  has  been  varied  linearly  between 
times  of  neutral  stability  (near  sunset  or  sunrise), 
(e)   N9  as  a  function  of  height  and  other  variables 

In  Eq .  (21),  it  was  shown  that, 

and  '..;e  now  have  convenient  working  equations  for  x  and  X, 
Eqs.  (34),  and  (39),  respectively.   Computations  of  N9  as  a 
function  of  elevation  can  now  be  made  provided  the  roughness 
parameter  z0  and  the  wind  profile  parameter  9   have  first  been 
obtained. 
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3.   Nature  of  the  data 

In  addition  to  testing  the  micrometeorological  theory 
forth  in  this  paper,  this  investigation  was  also  initiated  by  the 
availability  of  comprehensive  round-the-clock  data  of  the  highest 
caliber  to  date.   The  Great  Plains  Turbulence  Program   £7}  was 
a  research  project  of  the  U.  S.  Air  Force  Geophysics  Research 
Directorate  conducted  in  open  prairie  country  near  O'Neill,. 
Nebraska  between,  1  August,  and  10  September  1953. 

Fourteen  organizations  participated  in  the  project,  and 
their  contributions  included:   instrumentation,  data  collection 
and  evaluation,  development  of  turbulence  theory,  instrument 
description  and  evaluation,  site  preparation,  and  data  tabulation. 

The  Great  Plains  Turbulence  Program  was  in  progress  for  about 
six  weeks  during  which  time  seven  observation  periods  were  com- 
pletes.: .   The  criteria  for  an  acceptable  observation  period   ere 
that  the  daytime  periods  be  clear  to  partly  clou.^>,  with  clear 
skier,  at  night,  so  that  large  diurnal  variations  in  the  heat  budget 
terms,  and  particularly  in  the  Pvichardson  number  could  be  realize( 
Also,  it  was  preferred  that  the  mean  geos trophic  wind       vary 
from  one  observation  period  to  the  next.   Instrument  arrangement 
on  the  field  site  required  that  the  general  observation  periods 
be  characterized  by  relatively  constant  wind  direction  between 
3E  and  SW  £73  ■ 

Only  the  data  of  the  sixth  and  seventh  general  ol  ;ervation 
periods,  for  the  days,  31  August  through  1  September,  and  7-8 
tember  1953  CO    >  respectively,  will  be  described  and  evaluat 
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in  this  paper. 

Ma      iet.ts  Institute  of  Technology,  a  ad  Johns  H 

University  (hereafter  r<  Eerred  to  as  MIT,  and  JHU,  respectively), 
wind  and  temperature  profile  observations  at  standard  levels  were 
utilized,   The  MIT  data  collected  at  the  2,  4,  8,  and  16-meter 
levels  consisted  of  15 •minute  mean  wind  speed  (cm/sec),  and  13" 
minute  mean  air  temperatures  (°C).   JHU  data  was  collected  at 
,4,  ,3,    1.6,  and  3.2-meter  levels  and  consisted  of  hourly— mean 
wind  speed  (cm/sec),  and  average  temperatures  from  20  readings 
during  a  5-minute  interval.   Bi-hourly  observations  were  taken 
from  the  time  the  general  observation  period  began  with  all  time 
averag(   centered  at  observation  time. 

The  standard  cup  anemometers  used  by  MIT  were  of  the  conven- 
tional type  and  were  wind-tunnel  calibrated  prior  to  installation. 
Differences  in  calibration  were  minimized  by  field  matching  tests 
conducted  at  the  test  site.   Wind  measurements  were  averaged 
over  l:>-minute  intervals  and  recorded  to  the  nearest  cm/sec. 
JHU  obtained  a  number  of  three  cup  anemometers  of  1941-1942  model 
from  the  Soil  Conservation  Service,  and  by  modification  and  careful 
matching  produced  a  set  of  mechanical  wind  measuring  instrument:, 
with  compatability  of  the  order  of  0.5%,   Wind  measurements  were 
of  hourly  mean  wind  speed  recorded  to  the  nearest  cm/sec. 

MIT  temperatures  were  measured  with  identical* artificial ly 
ventilated  thermocouples,  similar  instrumentation  having  previously 
been  usee  at  MIT's  Round  Hill  Field  Station.   The  thermo junctions 
*  re  shielded  from  solar  and  ground  radiation  and  ventilated  at 

a  rate  of  5  m/sec.   Calibration  was  conducted  in  water  ath 
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matched  against  precision  mercury  thermometers.   Fifteen-minute 
mean  temperatures  at  the  standard  levels  used  in  this  paper  were 
recorded  to  the  nearest  Q.01°C. 

The  description  of  JHU  temperature  installation  is  similar 
to  that  of  MIT's,  the  main  differences  being  in  the  method  of 
shielding  and  the  lack  of  artificial  ventilation.   Calibration 
was  similar  to  HIT  procedure  previously  described,  Uiing  Beckmann 
differential  mercury  thermometers  between  the  two  water  baths, 
an  absolute  accuracy  of  0.05°C  was  indicated.   An  average  of  20 
temperature  readings  at  each  level  during  a  5— minute  period  was 
recorded  to  the  nearest  Q.Ol'C,   A  detailed  description  of  the 
instrumentation,  etc.,  is  given  in  Lettau   C73. 
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4 .   Computations 

The  solution  of  the  working  equations  for  ^  ,  zQ,    x,  and  X; 
equations  (8),  (11),  (34),  and  (39)  respectively,  was  accom- 
plished by  desk  calculator,  rather  than  by  slide  rule  or  electro- 
nic computer.   It  is  felt  that  use  of  the  desk  calculator  was 
justified  considering  the  relatively  small  quantity  of  data. 
Also,  greater  accuracy  than  that  attainable  by  slide  rule,  was 
assured,  to  the  limits  of  the  data  significance. 
(a)   The  roughness  parameter,  zQ 

Determination  of  z0  by  means  of  a  fitted  straight  line  on  a 
in  z   versus  wind  speed  graph  can  afford  only  a  rough  approxima- 
tion, at  best.   Accordingly,  the  working  equation, 


-c 

applicable  at  the  adiabatic  times  of  sunrise  and  sunset,  was  used, 
Fortunately,  the  0635  and  183 i>  observation  times  were  near  enough 
to  local  sunrise  and  sunset  at  this  time  of  the  year  (31  Aug.  - 
8  Sept.  1953)  to  satisfy  the  adiabatic  requirement. 

For  the  MIT  data  levels  of  2,  4,  8,  and  16  m,  Eq .  (11)  with 
Z4  =2^,  Z3  =  23,  Z2  =  2    ,    and  z^  =  2  meters  may  be  reduced  to 
the  form 


--  jmrca-5"MHy    (ho) 


The  roughness  parameter  for  the  JHU  data  levels  of  .4,  .8,  1.6, 
and  3.2  in,  is  obtained  by  reducing  Eq .  (11)  to, 

whe re  aga in ,  2  U=  U  H  +■  U4  -r  Ul+  U4  and ,  A  Q-  U<  -U-^  U^  -  U v 


:cl  by  reducing  hq .  (11;  to, 
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The  above  formulas  are  solved  very  simply  by  desk  calculator 
and  logarithmic  tables.   The  value  of  zQ   obtained  by  thi;;  method 
is  made  to  hold  for  all  computations  within  a  12  hour  period 
centered  on  the  time  for  which  it  is  determined. 
(b)   The  iterative  procedure  for  yi  =  1  -  6 

f£q .  (8)  required  only  that  the  roughness  parameter  be  known, 
in  addition  to  the  observed  data.   No  suitable  objective  method 
for  obtaining  an  initial  estimate  of  yj_  could  be  devised.   Accord- 
ingly, after  some  experience  in  the  use  of  iterative  steps,  it 
was  found  that  subjectively-determined  starting  values  for  the 
other  bi-hourly  times,  based  on  the  premise  that  ft  increases 
with  diurnal  heating  and  approaches  one  nearer  the  ground,  provided 
for  more  rapid  convergence  of  the  computation.   It  was  observed 
that  the  solution  converged  quite  nicely  for  all  but  near-neutral 
conditions,  where  considerable  difficulty  was  encountered. 
Accordingly,    :riteria  for  accuracy  of  the  solution  were  devised 
which  simplified  the  computations  a  great  deal.   That  is,  if 
ly.^1/  ,|0  ,  then  the  solution  could  be  halted  when  two  places 
of  decimals  had  been  fixed.   However,  for  l<,*J^.  .10  ,  and  if 
the  solution  was  oscillating  over  the  desired  root,  it  was 
required  that  three  places  of  decimals  be  unchanged  after  further 
iLerants  are  obtained.   If  the  solution  was  not  oscillating, 
but  approaching  the  root  from  above  or  below,  then  four  places 
of  decimals  were  determined  without  change  in  a  successive 
iterant.   Also,  it  'was  decided  to  disregard  the  solution  if 
lu.j|  <i •  005  ,  for  this  case  was  considered  identical  to  Q    =1. 
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In  addition,  certain  data-series  of  near  neutral  stability 
ro  excluded:   in  particular,  those  case.;  in  :  hich  It^-^M 
in  atiy  layer  was  less  than  0.05°C. 

(c)  Computation  of  x 

As  was  pointed  out  in  the  discussion  following  Eq,  (31) 
the  solution  for  x  requires  the  knowledge  of  ;:G,  a  ,  and  some 
representation  of  wa  under  non-neutral  conditions.   This  repre- 
sentation was  accomplished  in  Eq .  (34)  by  allowing  the  function 
AUc=(^4  "  Ma^J  +  \ty*,~  Uc'v,  )   lo  vary  linearly  between  the 
0635  and  1835  adiabatic  conditions.   The  \alue  of  x  obtained 
is  generally  negative  for  unstable  cases  and  positive  for  stable 
cast;, . 

The  roughness  parameter  has  already  been  determined  for  the 
solution  of  S  in  any  layer,  and  is  applied  in  the  solution  for 
x  in  the  same  manner  as  before.   The  value  of  -]_,2  appearing 
over  z0  iii  Eq .  (34)  makes  the  value. of  x  obtained  applicable  at 
I  .:■   geometric  mean  of  the  layer  concerned,  whereas  the   S  value 
was  assumed  to  be  representative  of  the  entire  layer  z\    to  zo. 

( d )  C  o mputationofX 

The  solution  of  Eq .  (39)  for  the  Richardson  number  for  any 
layer  is  quite  straightforward.   As  observation  shows,  it 
contain:,  only  constants  and  values  already  obtained,  or  from 
the  data.   The  value  for  g  was  taken  to  be  380  cm/ sec  .   To 
determine  the  value  of  02  -  ®i>    the  expression   <^z  -  ® t  '=t 
(Ti-T.)  t-   C2z-2,  )  &        was  used  where  £  ^  *C'a  ^tC:,; 
As  was  mentioned  earlier,  a  .criterion  for  determining  a  non- 
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neutral  layer  from  the  temperature  sounding  was  set  up,  namely, 
JO^-©,  J  s     .05°C.   In  the  contrary  case,  the  data  was  not  used 
for  computation  purposes.   The  value  o£  X  obtained  has  the  ;~.ame 
sign  as  the  temperature  lapse,  negative  for  unstable  conditions, 
and  positive  for  stable  conditions.   Since  z\.2   appears  in  the 
equation  for  X,  its  solution  is  also  applicable  at  the  geometric 
mean  of  the  layer.   The  working  equation  for  X  now  becomes, 

Y    -   ' &  *■** 3 ftT- "^  +  U~ 2' * ('%>■%*•*  ff  2.) 

The  solution  for  Nq  is  simply  the  value  of  x  from  Eq .  (34) 
divided  by  the  value  of  X  from  Eq .  (39).   The  value  for  Nq 
thus  obtained  is  always  positive,  and  applies  at  the  geometric 
mean  of  the  layer,  z±    to  :--2  • 
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5.   Discussion  of  results 

This  section  contains  a  summary  of  results  for  the  computa- 
tions of  6  ,  x,  X,  and  Nq,  and  compares  these  results  to  tho.se 
obtained  by  previous  investigators.   The  numerical  values  of  the 
above  parameters  are  shown  in  Tables  1,  2,  and  3  for  the  indicated 
times,  and  layers. 

Deacon  £2~)   quotes  values  of  the  wind  profile  parameter,  which 
he  obtained  from  Eq .  (5)  over  short  grass,  using  wind  speeds  for 
heights  between  the  0.5  and  h-n\  levels,  as  ranging  from  0.72  to 
1.13.   Inspection  of  Tables  1,  2,  and  3  shows  a  range  of  £  of  0.85 
to  1.14,  a  somewhat  lesser  range  than  that  obtained  by  Deacon  [_T\ 
from  the  series  of  data  he  utilized  at  Porton,  England. 

Priestley  [93  ,  shows  a  plot  of  log  Nq  versus  x,  for  mainly 
unstable  periods,  in  which  values  of  Rf  range  from  0.1  to  -0.5. 
His  scatter  diagram  is  shown  below  as  Fig.  1  of  this  paper. 

Based  on  this  data,  Priestley  obtained  an  average  coi rela- 
tion coefficient  between  Nq  and  Rfi  5  of  -0.81. 

1 


t 


N  =  V 


M 


Fig.  1 


4-0,1 


Dependence  of  K^/KM  on  Rf  at  1.5  meters 
(after  Priestley   l9j    ). 
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Table    1.      Computations    from   the   sixth  general   observation  period 
(MIT  data   only). 


u       ,  &  Time  of   observation 

2£    2£         31  August   1953 

T"  £e         0435  0835        1035  1235          1435          1635       2035          2235 

£      .9174  1.1274  1.0775  1.1436   1.0224      1.0145      .8692        .9027 

x      .5244  -.4987  -.3647  -.5385    -.2269                     1.1614 

X      .0799  -.0251  -.0860  -.0596   -.0888                       .0897 

Ne  6.5600  4.2420  9.0340  2.5560  12.9480 

f      .9922  1.1207  1.1093  1.1270   1.0714      1.0153                          .9823 

x  2.0605  -.4303  -.3891    -.3770      -.2949  .1629 


vO 


i 

CO 


co 

<*       X      .0267  -.0699  -.0738   -.0671      -.0341                          .0309 

Ne  2.2660  6.1540  5.2690   5.6270     8.6550                        5.2720 

S  1.0714  1.1173  1.0962    1.0475 

x  -.3226  -.4257  -.3900   -.2967 

™       X  -.0378  -.0745  -.0810   -.0625 

Ne  8.5420  5.7100  4.8160  4.7460 

1  September  1953 

0035  0235  0435  0635   0835 


\D 


i 

CO 


p 


9121  .8547           .9314        .9329      1.0674 

x      .6327  1.2431          .3809        .4000      -.3755 

X      .0916  .1007          .0804        .0622      -.0382 

NQ  6.9050  12.3450       4.7350     6.4280     9.8210 

(3     .9906  .9815 

x      .1044  .1837 

CO 

^   X   .0594  .0384 

N9  1.7560  4.7890 
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Table  2.   Computations  from  the  seventh  general  observation  period 
(MIT  data  only). 


Time 

Layer 

<? 

X 

X 

N9 

1235 

2-4 

1.0562 

-.9350 

-.0924 

11.118 

1435 

8-16 

1.1041 

-.5013 

-.0518 

9.669 

ii 

4-8 

1.0957 

-.9463 

-.0598 

15.819 

ii 

2-4 

1.0050 

-.9160 

-.0924 

9.909 

Table  3.   Computations  from  the  sixth  general  observation  period 
(JHU  data  only). 


Time 

Layer 

P 

X 

X 

N6 

(m) 

1035 

1.6-3.2 

1.0394 

-.2192 

-.0365 

5.999 

1235 

ii 

1.0657 

-.3085 

-.0493 

6.262 

1035 

.8-1.6 

1.0299 

-.1247 

-.0380 

3.280 

1435 

1.6-3.2 

1.0490 

-.2979 

-.0683 

4 .  365 

Table  4.   Computations  of  the  roughness  parameter 


Date  Time  zq   (™)  (MIT  data  only) 


31  August  0635  0.017960 

1835  0.012537 

1  September  0635  0.015483 

7  September  0635  0.174535 

1835  0.01368^ 


Date  Time  z.Q     tt»)  (JHU  data  only) 


31  August  0635  0.008249 

1835  0.007335 
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A  similar  comparison  made  from  the  results  of  this  investi- 
gation, shown  in  Figs.  2-6,  indicate  that  the  range  of  Nq  is 
several  times  as  great  as  that  shown  by  Priestley  L*9D, ^r  as 
calculated  by  Lettau  [jj   from  the  Great  Plains  data.   Likewise, 
values  of  x  computed  in  this  paper  ranged  up  to  -  1.0.   Exam- 
ination of  Priestley's  scatter  diagram,  Fig.  1,  shows  that  his 
range  of  x  <"0  is  much  smaller  (closer  to  neutral  conditions) 
than  the  points  for  which  x<0  in  Figs.  2-6.   There  is  also  an 
indication  in  Fig.  1  that  his  curve  is  making  a  bend  in  the 
vicinity  of  x=0,  similar  to  that  on  the  right  side  of  Figs.  2, 
3,  and  4. 

It  was  noticed  by  this  investigator  that  the  computation 
of  x  was  quite  sensitive  to  erroneous  values  of  the  roughness 
parameter.   In  the  case  of  very  rough  surfaces,  the  vertical 
wind  profile  may  be  represented  in  the  form 

-    jk  ^ 

where  d  is  called  the  datum-level  displacement  £3} .   Lettau  CZJ 
makes  use  of  the  zero— point  displacement  D=d  +  z0  in  computing 
values  of  zQ   at  the  JHU  site,  and  obtains  an  average  value  of 
d=9.9  cm.   This  means  that  the  logarithmic  profile  produces  a 
zero  value  of  wind  9.9  cm  above  the  zero  point  of  nominal  height, 
z0.   In  the  case  of  very  dense  vegetation,  d  +  z0  is  equal  to 
the  height  of  the  vegetation.   Lettau  [7]  reports  the  value  of 
d  he  obtains  to  be  on  the  order  of  the  grass  height,  10  cm. 
Lake  f5j  ,  and  Deacon  £2^   also  make  use  of  D  in  computing  values 
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of  z0.   Deacon  gives  values  of  z0  for  mown  grass~surfaces  of 
grass  length  4.5  cm,  ranging  from  zQ  =  2  A   cm  to  z0  ■-    \  J   era, 
depending  upon  the  strength  of  the  wind  at  the  2~m  level.   The 
results  of  this  investigation  gives  values  of  zQ  ranging  from 
0.73cm  to  l„80cm,  as  shown  in  Table  4. 

Lettau   ["7 3  calculated  values  of  X  using  an  equation  similar 
to  Eq»  (38)  at  the  2S4,  and  8-m  levels .   Values  of  X  obtained 
in  this  work  are  consistent  with  his  results  and  are  valid  at  the 
geometric  mean  of  the  layers  between  2,  4S  8S)  and  16  m„ 

Results  of  correlations  of  Ng  with  x,  X,  and  Q     ,  in  the 
entire  layer  2  m  to  16  ra  are  shown  in  Table  5. 

Table  5.   Correlations  with  Nq  with  x,  Xs  and  6  . 

Sample         Correlation  Coefficients         Significance 
size  Nq  vs  x         N0  vs  X        Ng  vs  ^       Level 

Stable   10   .965*         . 707*         .911*        .639 
Unstable  21   -.805*         -.130         -.244        =.440 


Note;   Significant  correlations  for  the  0.05  level  of  belief  are 
indicated  by  the  asterisk  (*). 

The  5%  significance  level  was  derived  from  the  test  men- 
tioned on  pp.  124  of   [41    All  of  the  correlations  obtained 
under  stable  conditions  were  positive  and  significant. 9  indicating 
that  N0  varies  with  stability  in  the  same  way  with  xs  X,  and   # 
However,  the  correlation  was  highest  with  xs  which  is  to  be 
expected  since  x  explicitly  contains  Ng  in  its  formula „   In  un= 
stable  layers ,  all  of  the  correlation  coefficients  were  negative 
although  only  one  of  these  was  significant  at  the  5%   level. 
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This  particular  correlation  verifies. in  a  general  way^that  of 
Priestley   ^93  ,also  taken  under  unstable  conditions „ 

Figo  7  shows  a  plot  of  Ng  versus  time  for  each  layer  in 
which  N©  was  determined.   For  the  period  31  August  -  1  September 
the  average  values  of  Uq   are  listed  below. 


Layer  (m)  Nq  ave 

8  -  16  7=5576 

4-8  4.9735 

2-4  5.9535 


Note  that  Priestley's  [9]    mean  value  of  Ng  at  105  m  would  be  near 
ly   one.   However j  he  has  values  of  Nq  as  high  as  4»0» 
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b.   Conclusions 

The  solutions  obtained  in  this  investigation  are  subject 
to  rather  stringent  assumptions  which  undoubtedly  account  for 
some  variance  of  the  results  from  those  of  previous  investigators 
It  has  been  shown  in  Section  5  that  values  of  z0«  ft  ,  and  X,  lie 
in  the  range  of  previous  investigations,  but  Ng  is  somewhat 
larger  than  comparable  values  previously  obtained* 

The  theory  presented  in  this  paper  is  developed  to  allow 
computations  of  the  various  parameters  for  non-adiabatic  condi- 
tions ,  since  the  'accuracy  criteria  of  Section  4  ruled  out  all 
computations  for  near-neutral  conditions =   By  Eq i,  (34 )s  errors 
in  ALi(X    which  tend  to  increase  x  tend  also  at  the  same  time  to 
decrease  X  [from  Eq ,  (39)3  ,    and  therefore  increase  Ngo 

However,  this  investigation  did  not  attempt  to  verify  the 
above  observation,  but  to  put  forth  a  method  of  computing  the 
pertinent  turbulence  parameters.   The  technique  of  representing 

ur,yji   UY  linear  variation  of  the  function  A  M.a  is  a  previously 
untested  one  and  may  be  the  contributing  factor  to  the  varia- 
tions of  Nq  thus  obtained. 

Since  these  computations  ignored  the  use  of  a  zero-point 
displacement  in  computing  values  of  the  roughness  parameter, 
no  conclusions  as  to  the  contribution  of  zQ  to  values  of  Ng  can 
be  made.   However,  reasonable  solutions  for  z0  were  obtained^ 

The  plot  of  Nq  versus  x  in  Figs,  2-6  shows  that  Ng 
increases,  both  with  degree  of  stability  and  with  degree  of 
instability,  and  is  near  one  for  near  neutral  conditions,,   On 
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the  basis  of  this  paper  it  seems  likely  Fig„  1.  would  also  show 
that  Ng  increases  with  stability  if  more  points  to  the  right  of 
the  origin  had  been  obtained„ 

The  correlation  of  Nfl  and  x  showed  a  higher  degree  of  sig- 
nificance than  that  of  Nq  versus  X.  This  was  to  be  expected  as 
N@  is  explicitly  contained  in  x„   D>ee  Eq.  (21)3  «   Priestley 

[9~]   also  reports  that  spuriously  low  correlations  of  N@  and  X 
can  be  expected  as  errors  in  measured  gradients  affect  Nq  and 
X  in  an  opposite  sense. 

Consistent  with  the  approach  of  this  investigation,  Section 
4  criteria  ruled  out  solutions  for  many  near-neutral  layers. 
This  investigation  revealed  that  the  seventh  observation  period 
(MIT  data),  and  JHU  data  for  the  sixth  observation  period  indicated 
near~neutral  conditions  much  of  the  time;  also  solutions  were 
unobtainable  in  the  first  five  observation  periods.   Therefore s 
it  can  be  concluded  from  this  investigation  that  only  the  sixth 
observation  period  of  the  Great  Plains  Project  contains  data 
giving  a  significant  range  in  stability,  for  the  purposes  of 
this  research. 

An  inspection  of  Fig.  7  shows  that  Nq  is  largest  in  the  highest 
layei   throughout  most  of  the  day.   Secondly,  N@  increases  with 
increasing  departure  from  adiabatic  conditions  in  either  sense. 
This  accounts  for  the  JHU  data  giving  fewer  results  at  its  lower 
observation  levels. 
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